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COLLIS 10NAL  DR  I  FT  WAVES  IN  A  PLASMA 


WITH  ELECTRON  TEMPERATURE  INHOMOGENEITY 

J.  F.  Drake  and  A.  B.  Hassam 
Department  of  Physics  and  Astronomy, 

University  of  Maryland, 

College  Park,  Maryland  20742 

ABSTRACT 

A  fluid  theory  of  collisional  electrostatic  drift  waves  in  a 
plasma  slab  with  magnetic  shear  is  presented.  Both  electron 
temperature  and  density  gradients  are  included.  The  equations  are 
solved  analytically  in  all  relevant  regions  of  the  parameter  space 
defined  by  the  magnetic  shear  strength  and  the  perpendicular 
wavelength  and  explicit  expressions  for  the  growth  rates  are  given. 
For  shear  strengths  appropriate  for  present-day  tokamak  discharges 
the  temperature  gradient  produces  potential  wells  which  localise  the 
mode  in  the  electron  resistive  region,  well  inside  the  ion  sound 
turning  points.  Mode  stability  arises  from  a  competition  between  th 
destabilizing  Influence  of  the  time  dependent  thermal  force  and  the 
stabilizing  influence  of  electron  energy  dissipation.  Convective 
energy  loss  is  not  Important  for  shear  parameters  of  present-day 


fusion  devices. 
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I.  INTRODUCTION 

Low  frequency  density  fluctuations  have  been  observed  in  a  wide  vari¬ 
ety  of  magnetic  confinement  systems  Including  tokaraaks,  stellarators,  and 
raultipoles . *  ^  Such  fluctuations  are  important  since  they  may  enhance  the 
crossfield  particle  or  energy  transport.  The  anomaly  in  the  crossfield 
thermal  transport  coefficient  has  been  most  extensively  documented  in 
tokamak  discharges  although  no  causal  relationship  with  the  density 

g 

fluctuations  has  been  established. 

It  is  widely  accepted  that  the  observed  density  fluctuations  arise 
from  drift  Instabilities  driven  by  the  expansion-free  energy  associated 
with  pressure  gradients  in  confined  plasmas.  Extensive  theoretical  Inves¬ 
tigations  of  drift  wave  stability  have  been  undertaken.  In  particular,  the 
influence  of  magnetic  shear  on  drift  waves  has  been  the  subject  of  substan¬ 
tial  controversy  for  nearly  a  decade.  Pearlsteln  and  Berk  found  that  shear 
caused  an  outward  convection  of  wave  energy  away  from  the  mode  rational 
surface  (£*B  =  0)  and  they  argued  that  drift  wave  stability  resulted  from  a 
competition  between  this  outward  convection  of  energy  and  the  release  of 

9 

expansion-free  energy  near  the  rational  surface.  Later  numerical  and 
analytic  calculations,  while  verifying  the  outgoing  wave  character  of  drift 
waves,  demonstrated  that  pressure  driven  drift  waves  are  always  stable  in  a 
sheared  magnetic  field  in  slab  geometry.*® 

The  majority  of  investigations  of  drift  wave  stability  have  allowed 
gradients  In  the  ambient  density  but  not  in  the  electron  temperature.  The 
dominant  motivation  for  the  study  of  drift  waves,  however,  is  the  experi¬ 
mentally  measured  anomaly  in  the  crossfield  thermal  conductivity  of 
electrons  in  tokamak  discharges;  In  these  discharges,  the  temperature 
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gradient  typically  exceed3  the  density  gradient.  Thus,  the  inclusion  of  a 
temperature  gradient  in  studies  of  drift  wave  stability  is  of  critical 

importance.  In  fact,  in  cecent  numerical  calculations,  the  temperature 

11  12 

gradient  was  found  to  alter  the  entire  structure  of  the  mode.  ’  It  was 
found  that  the  drift  wave  no  longer  had  the  character  of  an  outward  propa¬ 
gating  wave  but  remained  localized  in  the  vicinity  of  the  mode  rational 

Q 

surface,  well  within  the  "ion  sound  turning  points."  The  outward  convec¬ 
tion  of  the  wave  energy  therefore  does  not  have  a  stabilizing  influence  on 
the  mode  if  a  temperature  gradient  is  present.  This  localization  by  the 
temperature  gradient  appears  both  in  the  presence  or  absence  of  electron- 
ion  collisions,  v  .  In  the  collisionless  limit^  the  mode  remains  stable 
even  when  the  convective  energy  loss  can  be  neglected  because  >f  Landau 
absorption  of  wave  energy  by  electrons  in  the  vicinity  of  the  rational  sur¬ 
face.  In  the  collislonal  limit,  a)  <  v  ,  on  the  other  hand,  the  wave  can  be 

s 

destabilized  by  the  temperature  gradient  if  the  velocity  dependence  of 

12 

electron-ion  collisions  is  taken  into  account. 

In  the  present  paper  we  investigate  the  structure  of  collislonal  drift 
waves  in  a  sheared  magnetic  field  with  both  density  and  temperature 
gradients.  We  employ  a  fluid  treatment  which  complements  the  kinetic 
treatment  of  Ref.  12  but  which  allows  us  to  solve  for  the  frequencies  and 
growth  or  damping  rates  of  the  mode  analytically.  Our  conclusions  are  in 
qualitative  agreement  with  the  conclusions  of  Ref.  12.  The  temperature 
gradient  acts  both  to  localize  and  to  destabilize  the  drift  wave.  For 
parameters  typical  of  present-day  tokaraaks,  we  find  the  modes  to  be  strong¬ 
ly  localized  in  the  vicinity  of  the  mode  rational  surface,  well  inside  the 
ion  sound  surface.  As  the  shear  is  increased,  the  modes  become  leas  and 
less  localized  and,  eventually,  for  very  strong  shear  the  asymptotic 


behavLor  of  the  mode  far  from  the  rational  surface  assumes  the  outgoing- 
wave  structure  of  Pearlsteln  and  Berk. 

The  relative  stability  of  the  mode  results  from  a  competition  between 
the  destabilizing  Influence  of  the  "time  dependent  thermal  force"^  and 
local  resistive  energy  absorption  by  the  electrons.  Thus,  the  mechanism 

for  destabilization  of  the  drift  mode  Is  the  same  as  that  of  the  tempera- 

14 

ture  gradient  driven  ralcrotearlng  mode.  For  parameters  appropriate  for 
current  tokamak  discharges  (In  which  there  are  still  typically  many  modes 
with  to  <  )  these  drift  modes  may  be  unstable.  We  have,  however,  not 

included  finite  0  corrections  which  were  found  to  be  stabilizing  in  Ref. 

12. 

The  effect  of  temperature  gradients  on  collisional  drift  waves  was 
also  studied  by  Chen  et  al.1'5  Their  analysis  was  based  on  Braglnskii 's 
fluid  equations, ^  but  they  only  solved  the  equations  in  the  strong  shear 
limit.  Consistent  with  our  strong  shear  results,  they  found  that  the  tem¬ 
perature  gradient  does  not  localize  the  drift  mode.  However,  their 
solutions  differ  from  ours  in  this  limit  in  that  VTg  had  no  destabilizing 
influence  on  the  mode.  The  time  dependent  thermal  force,  which  is  the 

source  of  instability,  is  a  higher  order  correction  to  the  thermal  force 

1 3 

and  is  not  included  in  the  equations  of  Braglnskii;  hence,  this  force  did 
not  appear  in  their  equations. 

The  fluid  theory  presented  in  this  paper  employs  the  Braglnskii  trans¬ 
port  equations  modified  as  described  in  Ref.  13.  The  derivation  of  the 
basic  equations  describing  the  linear  drift  mode  are  presented  in  Section 
II.  In  Section  HI  the  structure  of  the  drift  wave  is  examined  as  a  func¬ 
tion  of  magnetic  shear,  as  represented  by  the  dimensionless  parameter 
S  =  (M/m)(L^,/Lg)  (u*-j./ve)  *  and  perpendicular  wavelength,  kpg,  where  M  and  m 
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are  the  Lon  and  electron  mass,  Lg  and  L^,  are  the  shear  and  temperature 

scale  lengths,  k  Is  the  perpendLcular  wavenumber,  pg  Is  the  Ion  Larmor 

radius  based  on  T  ,  and  w  =  kp  v  /L  is  the  diamagnetic  frequency.  The 
e’  *Ts8T 

structure  of  the  mode  falls  into  four  basic  categories  in  this  shear- 
wavelength  parameter  space  (see  Fig.  1) .  Eigenfunctions  and  eigenvalues  of 
the  drift  wave  are  calculated  analytically  in  each  of  these  categories.  A 
summary  of  our  results  appears  In  Section  IV  where  we  also  discuss  the 
relevance  of  the  results  to  present  magnetic  fusion  confinement  schemes. 


~V 
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-  1.61  T/ne^  ,  a  ■»  .71  , 

and  tildes  denote  perturbed  quantities.  The  perturbations  n,  T,  and  $  are 

2  2  2  2 

normalized  to  n,  T,  and  T/e,  respectively,  and  the  operator  V  -*-3  /3x  -k  • 

Tlie  quantity  a'  Is  a  numerical  coefficient  that  appears  as  factor  multiply- 

1 3 

ing  the  time  dependent  thermal  force.  An  accurate  and  reliable  value  of 

a'  has  not  yet  been  computed;  for  a  horentz  collision  model, ^  the  value  of 

1 8 

a'  is  105/16,  while  a  variational  calculation  including  electron-electron 

collisions  with  Z  «  ng/n^  »  1  yielded  the  value  .54.  Since  the  time 

dependent  thermal  force  is  the  mechanism  responsible  for  destabilizing  the 

present  mode  and  because  of  the  wide  disparity  in  the  values  of  a'  just 

quoted,  we  have  left  the  value  of  a'  unspecified. 

Equations  (l)-(5)  represent,  respectively,  the  equation  of  continuity, 

parallel  force  balance  for  electrons  with  the  term  proportional  to  w/ve 

being  the  time  dependent  thermal  force,  parallel  electron  heat  flow,  the 

parallel  component  of  the  curl  of  the  one  fluid  momentum  equation,  and  the 

parallel  one  fluid  force  balance  equation.  We  have  defined  the  parallel 

wavenumber  k„  ^  kx/L  . 

II  s 

It  will  be  3hown  later  that  except  for  the  case  of  very  strong  shear, 

(S  »  L  A.  ),  v ||  ^  can  be  neglected  in  Eqs.  (la)  and  (lc)  in  a  self- 

consistent  manner.  This  approximation  is  equivalent  to  discarding  the 

r 

"sound  term"  of  the  Pearlstein-Berk  equation.  Equations  (la)-(ld)  can  then 
be  reduced  to  the  eigenvalue  equation 

3^$/3s^  +  iV<j>  =  0  ,  (2) 


-■v 


where 
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SV 


b+ 


<532(L-*32) 

D(s) 


O  (u^T/u>)s 

d77) 


5  =  1  -  u>  /u>  , 

*n 


k  =  1.07  , 


0(3) 


1  +  (1+ic+y)s  +  xs 


S  =  (ui^/uJS  , 


b  H 


Y  5  (2/3) (1+a)  a  , 


2  2 

and  where  s  =*  [exp (i* /4 ) ] (x/A  ) ,  with  An  defined  by  k'  A  D  =  u> ,  and 

J  v  II  D  B 

D  =  T/(.51m  v  )  is  the  parallel  electron  diffusion  coefficient.  The 
ii  e 

independent  parameter  s  lies  approximately  along  the  line  exp  (-in /4)  with 
respect  to  the  real  x  axis.  We  will  solve  the  eigenvalue  equation  along 
this  rotated  axis  and  verify  that  the  asymptotic  solution  corresponds  to 
either  a  bounded  or  outgoing  wave  along  the  real  x  axis. 

Although  the  expression  for  the  potential  V  in  Eq .  (2)  appears  some¬ 
what  complicated,  its  general  structural  dependence  on  s  can  be  readily 
discerned.  In  fact,  for  practical  purposes,  the  true  potential  can  be 
replaced  by  the  model  potential  given  by 


V  <*  b 
m 


0 

~ 2 
(1+s  ) 


2 

s 


(3) 


The  actual  potential  given  in  Eq .  (2)  reduces  to  this  model  in  the  limit  x 
*  1,  Y  “  0.  In  the  absence  of  this  model  potential  has  two  turning 

points  (nearly  on  the  real  s  axis  for  5  <  0) ,  one  on  either  side  of  the 
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rational  surface  (s  =  0).  No  solution  localized  between  these  turning 

points  exists,  however,  because  the  potential  Is  evanescent  between  the 
2- 

turnlng  points  (V  $  -  0  along  the  real  axis).  The  only  possible  solution 
for  uj.rp  =  0  Is  a  solution  which  Is  bounded  between  this  Inner  turning  point 
and  the  outer  sound  turning  point.  In  the  presence  of  a  temperature 
gradient,  however,  the  foregoing  potential  has  four  turning  points,  two  on 
either  side  of  the  rational  surface.  A  localized  solution  between  each  of 
these  turning  points  Is  now  possible  If  6  has  the  appropriate  sign. 

Our  general  procedure  In  solving  Eq  .  (2)  Is  to  assume  that  ~ 

0(1),  and  to  examine  the  structure  of  the  potential  V  as  a  function  of  the 
parameters  3  and  b.  [In  most  magnetized  plasma  of  Interest,  and  are 
comparable  and,  In  keeping  this  ratio  fixed,  we  greatly  reduce  the  para¬ 
meter  space  and  hence  the  number  of  possible  parametric  regions.]  When  the 
shear  S  Is  very  weak,  we  find  that  the  two  turning  points  on  either  side  of 
the  rational  surface  are  nearly  coincident  (Region  I  of  Fig.  1).  The 
potential  well  then  becomes  approximately  parabolic  near  the  turning  points 
and  bounded  solutions  of  the  resulting  Weber  equation  can  be  obtained. 

These  solutions  correspond  to  the  bounded  solutions  obtained  analytically 
In  Ref.  12. 

/Vs  we  move  from  Region  I  >to  Region  II  of  Fig.  1,  the  two  turning 
points  begin  to  separate  and  the  inner  turning  point  moves  toward  s=0.  In 
this  region,  the  solution  of  the  Airy  equation  near  the  outer  turning  point 
Is  matched  to  the  solution  of  a  parabolic  cylinder  equation  in  the  vicinity 
of  s=0.  In  Region  III,  the  shear  S  is  large;  in  this  case,  the  potential 
Is  shallow  and  a  "constant  $  approximation"  can  be  made  to  solve  the 
equation.  This  constant  $  solution  near  s  -  1  is  matched  to  an  evanescent 
solution  for  s  »  l  which  Is  again  localized  well  within  the  Ion  sound 


turning  point.  FinaLly,  In  the  case  of  very  strong  shear  in  Region  IV  of 
Fig.  1,  the  "constant  solution  for  s  ~  I  is  joined  to  the  Pearlstein- 
Berk  outgoing  wave  solution.  In  this  region,  our  solution  is  consistent 
with  the  solutions  obtained  by  Chen  et  al.^  for  drift  waves  with  a 
temperature  gradient. 

The  detailed  solutions  of  Eq  .  (2)  are  described  in  the  next  section. 


[II.  ANALYTIC  SOLUTIONS 


As  previously  discussed,  Rq •  (?)  can  be  solved  in  the  four  regions  of 
parameter  space  defined  by  S  and  b  as  shown  in  Fig.  2.  We  now  examine  each 
of  the  regions  in  detail. 

A.  Region  I:  S  «  b3/2,l 

When  the  shear  is  very  weak  (i.e.,  as  S  -*■  0  )  the  mode  will  tend  to 
localize  Itself  in  a  region  where  the  right-hand  side  of  Eq.  (2)  Ls  small 
so  that  V(s)  *  S  1  remains  finite.  The  existence  of  such  a  localized  mode 
requires  a  second  order  turning  point  (3V/3s  *  0)  .  In  the  limit  S  -*■  0  , 
we  therefore  obtain  the  dispersion  relation  by  requiring  that  the  two 
turning  points,  given  by  solution  of  the  biquadratic  equation, 

x(b+6)s^  +  [  ( 1  -He  +y  ) b  +  6  -  o(u>  /co)]s2  +  b  -  0  ,  (4a) 

C  ’'It 


merge,  or 


[  ( I  -He  +Y )  b  +  6  -  o(wAT/u))]2  =•  4<  b  ( b+6  )  .  (4b) 

The  lowest  order  dispersion  relation  obtained  from  (4b)  is 

to  =  [co.  +  1.71  to  ]  (1  -  2.07  c,bl/2  +  .71ia'co_/v  )  b  «  1  , 

o  *n  *T  i  *T  e 

co  =  .28--  fl  +  .081a'  - — )  b  »  lf  n  »  1, 

o  b  ^  bv  J  e 

e 

1.71  n  1/2 
i  ^  1  +  1. 71  n  ‘ 


(5) 
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where  the  root  for  b  »  1  has  only  been  given  for  ne  >>  1  for  simplicity  of 
expression.  The  location  of  the  double  root  is  given  by 


2  <1/2  <  .  ..  . 

so  -  b  /c1  ,  b  «  1  , 

2  l  /? 

s^  =  {!/<)'  =  0.97  b  »  1  ,  ne  »  1  •  (6) 


Note  that  the  existence  of  a  double  root  requires  that  n  be  nonzero 

e 

(ng  >  1  for  b  »  1).  In  the  absence  of  shear,  there  are  no  stabilizing 
terms  in  the  expression  for  co  in  Eq.  (5)  and  so  the  time  dependent  thermal 
force  destabilizes  the  mode. 

When  S  is  finite,  the  roots  will  split  from  the  expression  given  in 
Eq.  (6)  by  an  amount  As  corresponding  to  a  shift  Aco  *  u)-wo  the  eigen¬ 
value.  By  expanding  V(s)  around  Sq,  we  obtain  the  eigenvalue  equation 

32$/3s2  +  [t — - ]  [  |  s-s  )2-(As)  2]^  =  0  ,  (7) 

S  sZ  D(s  )  0 

o  o 

whe  re 

{[(l+x-hr)b  +  6  -  0(uj  /id)]2  -  4xb(b-H5) }  ^2 

As  =  >TT7X -  "  •  (8) 

sks  (b+o) 
o 

The  bounded  solutions  to  this  equation  yield  expressions  for  the  shift  In 
the  eigenvalue  u>j. 


(tDjA oQ)  =  -  1.6  c^/2  exp(iit/4)  (S/b^2)  ^2  b  «  l  , 

1/2 

(u^/Wq)  =  -  1.6  exp(liv/4)S  '  b  »  1  ,  ne  »  1  ,  (9) 
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where,  again,  for  simplicity,  we  have  assumed  >>  1*  Since  Im  j)  ^  <  0, 
shear  lias  a  stabilizing  Influence  on  the  drift  wave  even  though  the  mode  is 
bounded.  This  damping  arises  from  electron  reslsltlve  dissipation  and  from 
parallel  heat  conduction.  The  Instability  drives  a  temperature  perturba¬ 
tion  [see  Eq .  ( 1 c) ]  which  Is  dissipated  by  parallel  diffusion 
along  B  and  Is  therefore  stabilizing.  The  overall  stability  of  the  mode 
results  from  the  stabilizing  Influence  of  this  local  dissipation  and  the 
destabilizing  Influence  of  the  time  dependent  thermal  force. 

The  localization  of  the  mode  Is  of  extent  As  around  sQ.  The  parabolic 
expansion  of  the  potential  well  around  therefore  requires  As  <<  sq  or 

S  «  b3/2  ,  b  «  1  , 

S  «  1  ,  b  »  1  .  (10) 

We  must  also  check,  of  course,  that  the  solution  obtained  from  Eq .  (7) 
Is  bounded  along  the  real  x  axis.  In  Fig.  2  we  show  the  anti-Stokes 
(solid)  and  Stokes  (dashed)  lines  for  the  potential  V  with  eigenvalues 
given  in  Eqs.  (5)  and  (7).  The  dashed  lines  are  constant  phase  lines  of 
the  elkonal  exp{ 1  /  ds'  [i  V(s')J  }  while  the  solid  lines  are  the  con¬ 
stant  amplitude  lines  of  this  elkonal.  The  solution  which  is  evanescent 
along  the  Stokes  line  emanating  from  the  turning  point  in  Fig.  2  is 
evanescent  throughout  region  B  to  the  right  of  the  anti-Stokes  lines  in 
this  figure.  Since  region  B  includes  the  real  x  axis,  our  solution  is  also 
bounded  along  this  axis. 

For  the  eigenvalues  given  in  Eq .  (5),  the  sound  turning  point  lies  on 

1  /2 

the  real  axis  at  |x|  *  x  (1+b)  with  x  -  w/k'c  ,  which  Is  well  outside 

s  sis 
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of  the  region  x  ~  A^  where  this  mode  Ls  localized.  The  neglect  of  the 
parallel  Lon  dynamics  Is  therefore  well  justified. 

B.  Region  II:  b3^2  <  S  <  1 

In  the  weak  shear  limit  previously  discussed,  the  mode  was  localized 
1/4 

near  ~  b  (for  b  «  1),  the  two  turning  points  being  separated  by  a 

distance  As  ~  (S/b^ ^ ) 3 ^ *  <<  sq.  As  the  shear  Is  Increased,  the  turning 

3  /2 

points  separate  until  S  ~  b  ,  As  ~  s^.  At  this  point  the  Inner  turning 
point  Sj  Is  near  the  origin  while  the  outer  turning  point  ls  around 

~  b3^  ~  <<  1.  In  Region  II  of  Fig.  1,  we  therefore  look  for  a 

solution  with  <<  s9  <<  1* 

When  h  <<  1,  the  positions  of  the  turning  points  s.  and  S2  can  be 
calculated  as 


s2  =  [a  (w^/(d)-A  ]  /5<  , 

s 2  =  b/ (5  k  s2)  .  (11) 

Defining  a  new  variable  p  i  3/82  and  assuming  s  <<  1,  we  can  rewrite  Eq .  (2) 

32l>/3p2  -  Q(p2-e2)(p2-lH  -  0  ,  (12) 


where 


2  2  2 
e  =  s^/32  «  l  > 

Q  =  -1  5k  s^/S  • 


(13) 


r 
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i 

i 


We  solve  this  equation  by  assuming  that  Q  Is  a  large  parameter,  and  match¬ 
ing  the  Airy  solution  of  Eq .  (12)  near  p  =  1  with  the  solution  of  the 

2  2 

parabolic  cylinder  equation  In  the  region  p  “  e  a  0  through  a  WKB 
region  e  <  p  <  1.  The  solution  for  $  near  p  “  1  which  is  decaying  for  p  > 
1  Is  given  by  $  =  A1 ( (2Q) ^ ^(p-1 ) ] .  ^  The  WKB  solution  which  matches  the 
p<l  behavior  of  this  function  Is  given  by 

♦  -  Bin  [Q1/2(l-p2)3/2/3  +  it /4]  .  (14) 


Near  the  origin  this  solution  breaks  down  but  Eq.  (12)  simplifies  to 


3  24>/3p2  +  Q(p2-e2H  -  0  , 


(15) 


which  Is  the  parabolLe  cylinder  equation.  We  are  actually  only  Interested 

2  1/2 

In  solutions  for  which  £  Q  <<  1,  i.e.,  the  Inner  turning  point  can  be 

2  1/2 

entirely  neglected.  That  £  Q  Is  actually  small  shall  be  demonstrated _a 

posteriori .  Equation  (15)  now  becomes  a  Bessel  equation  with  the 
19 

solution 


1/2  r 
P  c 


rnl/2  2. 

1/4  [Q  P  /2]  , 


(16) 


where  Cy  Is  Yy  for  4>  even  and  Jy  for  <f>  odd.  The  dispersion  relation  la 
then  easily  obtained  by  matching  the  asymptotic  solution  for  In  Eq.  (16) 
with  the  small  argument  limit  of  $  in  Eq.  (14).  We  find 


,1/2 


(n  +  1/4)3* /2  , 


(17) 
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where  u  Is  a  positive  Integer.  For  n  >  0,  the  solution  of  this  dispersion 
relation  yields  a  Stokes  structure  which  Is  similar  to  that  shown  In  Fig.  2 
except  that  the  Inner  turning  points  are  close  to  the  origin.  For  n  <  0, 
the  Stokes  structme  Is  entirely  different  and  our  evanescent  solution  for 
p  »  1  does  not  correpond  to  an  evanescent  solution  along  the  real  axis. 

The  dispersion  relation  In  F.q .  (17)  can  be  solved  for  the  eigenvalue,  viz., 

2  2  1/1 
w  =  +  1 . 7  1.jo*t  +  (1 .21/c^a '(ojJkT/ve)  -  c2  exp(lir/6)  S 

c  2  =■  1.81  (1  +  l.71n  )  (n  +  1/4) 2/3  cj/3  .  (18) 

In  the  absence  of  the  destabilizing  Influence  of  the  time  dependent  thermal 

1  /3 

force,  the  mode  Is  damped  at  a  rate  which  scales  as  S'.  The  location  of 
the  turning  points  can  be  calculated  with  the  expression  for  w  given  in  Eq. 
(18).  We  find  s2  ~  S1  ^ ,  while  Sj  ~  b^^/S3^.  The  condition  Sj  <<  S2 
then  requires 

3/2 

S  »  bJ/Z  ,  (19) 

2  1/2 

which  also  guarantees  that  £  Q  «  1.  The  assumption  s2  «  1  along  with 
the  requirement  Iiwa)  <<  Reu>  requires 

S  «  1  .  (20) 

Region  til:  1  <  S  «  Lg/LT 


i 


When  S  is  very  large,  the  potential  V  ~  S  ^  is  weak  and  the  drift  wave 
is  no  longer  localized  in  the  region  s  ~  l,  i.e.,  the  shear  prevents  strong 
localization  of  the  mode.  However,  in  Region  ill  we  can  still  obtain  solu¬ 
tions  which  are  localized  within  the  ton  sound  turning  point.  Equation  (2) 

1/2 

is  solved  in  the  two  regions  s  «  S  and  s  »  S  and  the  solutions  are 

1/2 

matched  in  the  overlap  region  S  <  s  <  S.  When  s  <<  S,  Eq.  (2)  can  be 
approximated  as 

3  2 / 3  s ^  ~t— - —  {  [1+(ic+y)s2]6  +  a(w^  /u>)s2}  $  ■  0  ,  (21) 

S  D(s) 

_2 

where  we  have  neglected  (b+6)/S  ~  S  (to  be  checked  self-consistently) . 
Since  S  is  large,  this  equation  can  be  solved  iteratively.  To  lowest 
order,  <fr  -  const.  Corrections  are  obtained  by  Integrating  Eq.  (21) , 

s  s '  || 

<|>(s)/4>(o)  =  1+1/  ds  '  /  {  [1+(k+y )  s" “ ]6  +  o(u>*  /w)  s"2} 

S  o  o  UKS  '  l 

-  1  +  i(s/S)[6I1  +  a(w*T/a))I2]  ,  (22) 


wiie  re 


1 1  =  f  ds  [l+(<+y)s2 ]/D(s)  =  5.0  , 
o 

00 

1?  I  {  ds  s2/D( s)  =  1.2  , 
o 


1/2 


and  the  second  step  in  (22)  is  valid  for  l  «  s  «  S.  In  the  region  S 
«  s,  Eq.  (21)  simplifies  to 


18 


3  '  4>/3  s  +•  1  { ( b+6 )  / S  ] 0  =  9  , 


(23) 


with  the  bounded  solution 


$(s)  ~  exp  {  [  ( b-HS )  /  S  ] 1  ^  2  exp  (3ni/4)}  . 


(24) 


Matching  the  two  solutions  In  Eqs.  (22)  and  (24),  we  find  the  dispersion 
relation 


(b-h$)1/2  =  e-1"74  S~172  [6^+  o(w*T/w)I2]  (25) 

and  the  corresponding  eigenvalue 

w  =  W(1  “  1  C3/S)  * 

c3  5  (4-4/neHne  '  2,4  l+b^2  *  (26) 

which  has  a  real  frequency  which  is  characteristic  of  the  drift  wave  with¬ 
out  np*  Note  also  that  this  mode  Is  always  stable.  In  addition  to  the 
requirement 

S  »  1  ,  (27) 

which  was  Invoked  to  justLfy  the  iteration  procedure,  we  require 

1/2 

Re(b+6)  >  0  so  that  in  Eq.  (24)  Is  bounded.  With  the  solution  given 

in  Eqs.  (25)  and  (26) ,  this  requires  >  b/(l+b)  so  that  this  solution  is 
only  possible  when  is  nonzero.  Finally,  we  require  that  the  mode  be 
localized  with  the  ion  sound  turning  point.  The  scale  length  of  the  mode 
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solution  as  s  -*■  0  is 


r 


l  3$ 


(32) 


Equating  this  expression  to  the  Logarithmic  derivative  of  the  solution  for 
the  inner  region  given  in  Eq.  (23),  we  obtain  the  dispersion  relation 


1/2 


61. 


+  a(u>*T/a))  I_]  . 


(33) 


When  S  «  (L,,/!^)  ,  the  solution  to  this  dispersion  relation  corresponds  to 

a  »  1  and  simply  reduces  to  Eq.  (26).  The  ion  sound  term  is  not  important 

in  this  limit.  When  this  Inequality  is  reversed  (Region  IV  of  Fig.  1), 

1/2 

r( 1/4  +  a/2)  ~  (S  L  /L  )  »  1  and  the  solution  is  given  by  a  *  -  1/2  or 

n  s 


(0 


_  r ! 

1+b  1 


4c,  1/2  L  1/2 

t-2-)  (t£)  l  • 

s 


(34) 


The  damping  term  in  Eq.  (34)  corresponds  to  the  usual  convective  energy 
loss  Ln  a  sheared  magnetic  field.  The  temperature  gradient,  on  the  other 
hand,  simply  causes  a  shift  in  the  real  frequency  of  the  mode  in  this 
strong  shear  regime.  The  destabilizing  influence  of  the  time  dependent 
thermal  force,  though  present,  is  small  In  this  limit.  The  dispersion 


relation  in  Eq.  (34)  corresponds  to  that  given  in  Ref.  15. 
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fit.  SUMMARY  AND  DISCUSSION 

Che  colLisional  drift  wave  in  the  presence  of  7T  falls  into  four 

r  e 

distinct  regimes  as  a  function  of  magnetic  shear  and  the  perpendicular 

wavelength.  These  four  regions  are  shown  in  the  S-b  phase  space  in  Fig.  1, 

2  2  2  2 
where  S  =  (M/m)  (L~/L  Is  the  dimensionless  shear  and  b  »  k,p  Is  the 

T  s  *1  e  is 

perpendicular  wavelength.  In  Region  I,  the  shear  Is  weak  and  the  modes  are 

strongly  localized  on  either  side  of  the  rational  surface  at 
,  2L  1  /2 

x  ~  +  A^  =  +(w/k'  D  n )  .  In  Region  II,  the  shear  Is  somewhat  stronger  and 

the  mode  penetrates  to  the  rational  surface  but  Is  still  localized  inside 

the  distance  |x|  <  A^.  In  Region  III,  the  shear  is  strong  enough  to 

prevent  strong  localization  of  the  mode  within  the  distance  |x|  <  A^,  but 

is  sufficiently  weak  that  the  mode  is  still  localized  well  within  the  ion 

-1  /2 

sound  turning  points  (|x|  ~  p  (L  /L  ) (l+b)  ].  Finally,  when  the  shear 

v  s  s  n 

i3  very  strong,  the  mode  assumes  the  usual  outward  propagating  structure 
characteristic  of  the  drift  wave  with  VT^  =  0.  In  Regions  I  and  II,  the 

real  frequency  of  the  mode  scales  roughly  as  (w  +■  w  )/(l+b)  and  is 

*T  *n 

determined  by  the  electron  dynamics  In  the  region  x  ~  A^.  In  Regions  III 

and  IV,  the  frequency  is  m  /(1+b)  and  is  given  by  the  plasma  dynamics 

ouslde  the  electron  dissipative  region  |x|  >>  Ag. 

The  relative  stability  of  the  collisional  drift  wave  results  from  a 

competition  between  the  destabilizing  influence  of  the  "time  dependent 

thermal  force"  and  the  stabilizing  influences  of  electron  dissipation  and 

convective  energy  loss.  In  the  limit  S+0,  the  mode  is  unstable  with  a 
2 

growth  rate  y  ~  u)._/v  .  As  the  shear  increases,  electron  collisional  dis- 

sipation  and  parallel  thermal  conduction  have  a  stabilizing  effect  (scaling 
1/2  1/3 

as  S  in  Region  I  and  S  in  Region  II)  and  the  growth  rate  of  the  mode 


decreases.  In  Region  III,  the  drift  mode  is  always  stable  with  electron 
dissipative  effects  completely  dominating  the  destabilizing  thermoelectric 
terms.  Finally,  in  Region  IV,  the  convective  energy  loss  becomes  important 
and  the  mode  is  again  stable. 

Since  the  eigenf requencies  of  the  drift  mode  have  been  calculated  for 

all  values  of  shear  and  wavelength  (under  the  restriction  u>/v  .  <<  1),  the 

ei 

stability  of  a  given  mode  for  a  given  set  of  parameters  can  be  calculated. 

It  would,  however,  be  useful  to  obtain  a  shear  stability  criterion  which 

guarantees  stability  for  all  modes  k^  for  a  given  shear.  To  do  this,  we 

2  2 

note  that  both  b  =*  k^pg  and  S  «  *  k^  are  functions  of  the  perpendicular 

wavelength.  Thus,  by  increasing  k^  for  a  fixed  shear  and  collisionality , 

we  may  trace  out  a  curve  in  the  S-b  phase  space.  The  two  classes  of 

possible  curves  are  shown  by  the  dashed  lines  A  and  B  in  Fig.  1  for 

1/2  2  > 

(M/ra)  ve/Lgv^  <  1,  respectively.  For  case  A,  the  only  possible 

unstable  modes  are  those  Ln  Region  II.  The  conditions  y  >  0  with 

S  <  1  and  w/v  <  1  yield  the  constraint 
el 


V. 


M 

(_ 

m 


L2  -1 

.  2j 


>  w/v 


which  can  only  be  satisfied  for  Lg  /L^, 
require 


M  LZ 

.i  >  t:  ri>  • 

u 

s 


>  (M/ra) 


i/2 


(35) 


Thus,  for  stability  we 


Lg/LT  <  (M/m)1/2. 


(36) 


This  stability  threshold  Is  rather  sensitive  to  the  value  of  a'  and 
therefore  should  only  be  considered  approximate.  In  Case  B,  the  modes  can 
be  unstable  In  either  Region  I  or  Region  II.  The  modes  have  larger  growth 


rate  In  Region  1  with  b  '>  1,  so  that  If  tile  drift  wave  is  stable  here  It 

wlLL  also  be  stable  In  the  other  regions  of  Fig.  1.  In  Region  I  with  b  > 

1,  the  condition  r  >  0  with  w/v  <  1  and  S  <  1  again  yields  the  constraint 

e 

on  <jj  /v  ^  given  In  Kq .  (I1));  so  all  modes  will  again  be  9table  when  the 
inequality  In  Eq •  (lb)  Is  satisfied. 

Relevance  to  T  »kaiaak _ P 1  srim rges 

We  now  t  iro  our  attention  to  tlie  Importance  of  the  VT  -modified 

e 

drift  w  i  ,v>s  to  t  >kanak  discharges.  Consider  modes  with  k^pg  -  1;  then,  the 

1  /2 

1  aequo  t  I  ty  -  '  '  k.o^c,/!.  <  requires  Lc/qR  <  (M/m)  (a/qR),  which 

Is  df ( ■  Suit  t  >  sat!sf-  : .  many  tokamak  discharges.  (L.^  is  the  collisional 
ne  in  ‘ ■  'it':.  '  !hus,  i  r  i  ft.  wave  specialists  have  tended  to  regard 

t  'Vatu'- s  is  be  l  ny  ••  >1  1  is  1  on  less  or  nearly  collistonless,  and  Instabilities 
such  is  r  ie  dissipative  trapped  •  •lectrou  Instability  or  the  collislonless 
'73  resonance-  !  r  1  .vn  irlft  wave  are  considered  to  be  more  Important  titan  the 
coll  isional  drift  wave.  On  the  other  hand,  there  are  still  large  numbers 
of  collisional  modes  present  even  In  the  highest  temperature  discharges. 

Tlie  condition  *)#  <  v  for  all  modes  in  a  tokamak  requires  the  Inequality 

1  /2 

C^/qR  >  (M/m)  (a/p ^ ) (3/qR ) ,  which  Is  far  from  being  satisfied  in  present 
machines.  There  Is,  at  present,  no  compelling  theoretical  evidence  for 
believing  that  collisional  drift  waves  are  more  Important  than  collislon¬ 
less  modes  or  vice  versa. 

Density  fluctuations  have  been  experimentally  measured  In  a  variety  of 
tokamak  discharges  by  a  number  of  different  techniques  including  microwave, 
COr,,  and  far  Infra-red  scattering,  and  probes.  The  spectra  are  character¬ 
istically  peaked  at  the  longest  measurable  wavelengths  (lowest  frequency) 
while  the  fluctuations  are  peaked  at  a  minor  radius  rg  ~  a/2.  In  Table  I, 
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we  summarize  the  d  it  i  if  some  published  measurements  of  density 
fluctuations  in  tokamaks.  The  first  five  columns  list,  respectively,  the 
tokamak  on  which  the  measurements  were  made,  the  position  (minor  radius)  of 
the  fluctuation  measurements ,  the  local  temperature,  the  local  density,  and 
the  local  Zgjf,  as  reported  by  the  authors.  In  instances  where  only  peak 
(r=0)  expressions  for  Te  and  ng  were  reported,  we  approximated  T^(r^)  - 
T  (0)/2  and  n  (r  )  ~  .75  n  (0).  In  the  sixth  column,  we  have  tabulated 

6  Q  S  (2 

the  90°  electron-ion  scattering  collision  frequency,  v  ,  as  computed  from 

the  local  plnsma  parameters.  Finally  in  the  last  column  is  the  spectral 

half  width,  w .  =  2irf  ,  of  the  observed  fluctuations,  where  fKm  is  the 
’  nm  hm  nra 

frequency  at  which  S (ui )  falls  approximately  to  half  of  its  peak  value. 

Note  that  u),  ~  v  for  all  the  reported  measurements.  It  should,  however, 
hm  <  ei 

be  emphasized  that  the  PLT  data  reported  in  Ref.  4  was  taken  before  neutral 
beam  heating  was  undertaken  and  that  this  data  may  therefore  not  be  repre¬ 
sentative  of  present  PLT  operation.  Nevertheless,  Table  I  indicates  that 
collisional  modes  are  quite  likely  to  be  more  important  in  present-day 
tokamaks  than  has  been  previously  believed. 

The  observation  of  collisional  modes  in  tokamaks  does  not  guarantee 
that  the  fluctuations  are  destabilized  by  the  mechanism  reported  in  this 
paper.  The  shear  stabilization  condition  reported  in  this  paper,  [Eq. 
(36)],  is  not  well  satisfied  in  present  tokamak  discharges  so  that  insta¬ 
bility  is  probable  within  the  limitations  of  the  present  model.  However, 

finite  3  effects  were  found  to  have  a  stabilizing  effect  on  this  mode  in 

1  2 

previous  numerical  computations.  The  effect  of  finite  0  on  the  present 
fluid  description  of  this  instability  is  not  yet  determined. 

Perhaps  the  most  important  conclusion  of  this  work  is  that  the  temper¬ 
ature  gradient  strongly  modifies  the  structure  of  collisional  drift  modes. 


-4 
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Colllslonal  modes  Ln  present  tokamak  discharges  would  fall  along  curve  A  In 
Regions  II  or  III  of  Fig.  1  so  that  the  modes  would  be  localised  and  the 
convective  loss  of  wave  energy  not  important.  This  conclusion  may  elimi¬ 
nate  the  original  motivation  for  studying  the  effect  of  toroidal  coupling 
on  colllslonal  drift  waves,  i.e.,  to  nullify  the  convective  loss  of  wave 
energy  by  coupling  adjacent  pololdal  modes.  However,  the  existence  of  a 
new  unstable  branch  of  the  colllslonal  drift  wave  in  a  4  0  plasma, 
analogous  to  the  "toroidicity-induced"  branch  for  '7T^  «  0  plasmas,  J  can 
not  be  discounted. 


i 

•i 


I 
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FICURK  CAPTIONS 

Fig.  1.  The  four  regimes  of  the  VT^-raodif led  drift  wave  as  shown  in  the 

2  2  2  2 

S-b  phase  space,  where  S  =  (Ml.  /mL  )(w  /v  )  and  b  =  k  p 
1  n  s  *T  e  L  s 

represent  the  magnetic  shear  strength  and  perpendicular 
wavelength,  respectively. 

Fig.  2.  Anti-Stokes  (solid)  and  Stokes  (dashed)  structure  for  the 

VT  -modified  drift  wave.  The  Ion  sound  turning  points  fall 
e 

outside  the  domain  of  the  plot. 
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